function triangle_grid(JM, JXY)
E = length(JM); % 读取节点信息矩阵的行数, 即单元总数
N = length(JXY); % 读取节点坐标矩阵的行数, 即节点总数
x = JXY(:,1); % 建立总的 x 向量
y = JXY(:,2); % 建立总的 y 向量
ele_num_print = 1; % 是否显示单元编号
pnt_num_print = 1; % 是否显示节点编号
axis equal;
hold on;
del_x = 0.001;
del_y = 0; % 单元号放置的位置调整变量
for k = 1:E
    for l = 1:3
        p = JM(k,l);
        xx(l) = x(p);
        yy(l) = y(p); % 读取一个单元中的对应点的坐标
    end
    xx(4) = xx(1);
    yy(4) = yy(1);
    plot(xx, yy);
    x_cen = sum(xx(1:3))/3;
    y_cen = sum(yy(1:3))/3;
    if (ele_num_print == 1) % 标记每个单元的单元编号
        text(x_cen-del_x, y_cen-del_y, ['(', int2str(k), ')']);
    end
end
if (pnt_num_print == 1) % 打印节点编号
    for n = 1:N
        text(x(n), y(n), int2str(n));
    end
end
axis off

end